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ABSTRACT 

HR1217 is one of the best-studied rapidly oscillating Ap (roAp) stars, with eight known 
oscillation modes that are distorted by a strong, global magnetic field. We have reanalysed 
the multisite observations of HR 1217 taken in 1986 and 2000. We determined a weighting 
scheme for the 1986 and 2000 data to minimize the noise level. A wavelet analysis of the data 
has found that the modulation of the amplitude due to rotation for all frequencies is, in gen- 
eral, consistent with the expected modulation for modified £ = 1, 2 or 3 modes. Unexpected 
variations in the rotational modulation are also seen, with variations in the modulation profile, 
time of maximal pulsation, and pulsational energy in each mode. Interestingly, these changes 
take place on a short timescale, of the order of days. We consider potential explanations for 
these behaviours. 
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1 INTRODUCTION 

Many stars are observed to oscillate with periods that depend on 
their internal structure. Characterising oscillations therefore al- 
lows stellar parameters to be constrained. Rapidly oscillating Ap 
(roAp) stars are particularly interesting because their strong mag- 
netic fields, of the order of a kiloga uss (up to 25 kG), can s ignif- 
icantly distort the oscillation modes jCunha & Goughll2000l) . Fur- 
thermore, the observed longitudinal magnetic field strength and the 
pulsation amplitude are modulated by the rotation of the star, al- 
lowing the mode geometry of ro Ap stars to be studied in a way not 
possible for other pulsating stars (lKurtJl982l) . 

The roAp stars pulsate in high radial overtone modes simi- 
lar to those seen in solar-like oscillators. With the radial overtone, 
n, much greater than the spherical degree, £, the mode frequen- 
cies are roughly equally spaced by the so-called 'large separation'. 
This asymptotic frequency spacing, Au, is a measure of the sound 
crossing time of the star, which in turn is determined by the star's 
mean density jGabriel et al.lll985l) . Assuming adiabatic pulsations 
in spherically symme tric stars, the pulsation frequencies are ap- 
proximately given by ( iTassoul|[l98bi : ,Gough. 1 98dh 



one of the first discovered roAp stars dKurt j| 1 982h . It was the sub- 
ject of an extensive photometric global campa ign in 1986 , whic h 
detected six principal pulsation frequencies dKurtz et al.l Il989h . 
Five of these were alternatively separated by 33.4 and 34.5 /iHz. 

Without a precise identification of the degree, £, of the pulsa- 
tion modes, the large separation is uncertain by a factor of two. If 
modes of alternating even and odd degree are present they will be 
separated by about Ai'/2 in frequency, but if only modes of the 
same degree and consecutive n are present, they will be separated 
by about Ai/. It was inconclusive from the 1986 observations as to 
whether the large separation was 34 or 68 /iHz. While a relation has 
been observed between Av and the frequency of maximum power, 
I'niax, for solar-like osci llations that can resolve this type of ambi- 
guity jStello et alj2009l) . such a relation has not yet been shown to 
exist for roAp stars. 

The ambiguity can also be resolved by a precise de- 
termination of the luminosity. A large separation of 34 /iHz 
would imply that HR 1 217 is well above the Main Sequence 



Heller & Kawaled Il988h . The Hipparcos parallax measurement 



Perrvman et al.l 119971) showed that HR 1217 is close to the 



Au{n + 1/2 + e). 



(1) 



HRI217 (DOEri, HD 24712) is one of the best-studied and 
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Main Sequence, and hence the lar ge separation must be 68 /iHz 
jMatthews, Kurtz & Martinet 1 199 9). In principle the ambiguity 
could also be resolved by observations with frequency resolution 
high enough to separate the modes with 1 = and £ = 2 and 
modes with £ — 1 and £ — 3, for then the different 'small separa- 
tions' i^n,e — u„-i,e+2 would reveal the values of £, as they did for 
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the first identification of the low-degree modes detected in whole- 
disc observations of the Sun. 

A sixth frequency in the 1986 data at ~ 2.806 mHz was 
separated by 50 /^Hz from the fifth which, at a separation of 
|Ai/, did not fit the asymptotic relation. However, that rela- 
tion is valid for linear adiabatic pulsations in spherically sym- 
metric stars, whereas strong magnetic fields and chemical inho- 
mogeneities in roAp stars significantly break the spherical sym- 
me try. The effects of che mical inhomogeneities were examined 
by iBalmforth et"ai] ( 1200 ll) . but of greater interest to us here is 
the effect of the magnetic field. To estimate this a singular per- 
turbation approach has been app lied to determine the modifi- 



cation to the frequency spectra ( Dziembo wski & Goodd [l996l; 
Bigot et al] l200d: ICunha & Goughl I200Q: iBigot & Dziembowskil 



20021 : ISaio& Gautschvll2004h . ICunha & GougH ( l2000h found that 
at some frequencies the magneto-acoustical coupling is partic- 
ularly strong, resulting in a loss of a mplitud e from the mode 
and a frequency shift of several /^Hz. ICunhj ( 1200 ll) suggested 
such a mechanism may be responsible for the unusual separa- 
tion of the sixth frequency of the 1986 data, and predicted an 
additional mode 34 ^Hz abov e the fifth frequenc y that had been 
missed by the observations of iKurtz et al. I( fl98^ . More recently, 
ISaio, Rvabchikova & Sachkovl i 2010l) suggested that this mode 
might correspond to a modified I = ?> mode. 

HR 1217 was again observed photometrically in November 
and December 2000 as a part of the Whole Earth Telescope (WET) 
extended coverage campaign (Xcov20), with telescopes from eight 
observatories with apertures ranging from 0.6 to 2.1m, in one of 
the highe st-precisio n ground-based photometric studies ever un- 
dertaken (Kurtz et al.;;2005j ). In a preliminary analysis of the data, 
the h ypothesized 'missing' mode was indeed found dKurtz et^ 
l2002h . with amplitude lower than those of the six modes that had 
been found already. From a more detailed anal ysis, two closely 
spaced modes at this frequency dKurtz et a l.'2005') were found. This 
finding appeared to support the theory of Cunha & Gough (200ol). 
Moreover, the greater damping expected from the theor y might pro- 
vide th e explanation for the lower amplitude. However. IShibahashil 
( l2003h expressed the concern that the magneto-acoustic coupling 
could result in damping so strong that the mode would probably 
not be excited, rendering the reason for the existence of the modes 
an open question. However, no calculation of the excitation of the 
modes has yet been carried out to justify or allay that concern, not 
even by a generalization of the simple, highly idealized, procedure 
adopted by Balmforth et al. (2001). 

It must be asked why the 'missing' modes were missing in 
the first place, given that they were found in the 2000 data with an 
amplit ude greater than the noise level of the 1986 data. lKurtz et al.l 
( l2005h noted that the amplitudes of the other modes had also varied 
between the two datasets by several tenths of a mmag, which is 
greater than the amplitude of the 'missing' modes in the 2000 data. 
Such a change in amplitude could well provide the explanation for 
why the 'missing' modes were originally absent. 

An interesting observation made by Kurtz et al.l ( |2005|) was 
that there was little difference in the sum of power between the two 
data sets, indicating that, since the observed modes all have about 
the same inertia, the total pulsation energy may have been con- 
served, suggesting that nonlinear interactions between modes trans- 
fer energy on a timescale shorter than the intrinsic growth times. 
Other roAp stars that have multiple pulsation modes show signifi- 
cant amplitude variati ons on timescales of a few days, H P 60435 
being one example dMatthews. Wehlau & Kurt3 Il987l ). On the 
other hand, pulsation amplitudes of singly periodic roAp stars, such 



as HR3831 dKurtz et^l I997h . appear to be very stable. However, 
some roAp stars appear not to exhibit conservation of total pulsa- 
tion energy. For example, power is clearly not being conserved in 
HD 217522, in w hich a new frequenc y appeared between data sets 
in 1982 and 1989 dKreidl et alj|l99lh . 

A unique aspect of the oscillations of roAp stars is that the pul- 
sation amplitude is modulated with the rotation of the star, as is the 
apparent magnetic field strength. The modulation is well explained 
by the oblique pulsator model, in which the pulsations are axisym- 
metric oscillations with the sa me axis as the magnetic field and in- 
clined to to the rotational axis dKurtzlll983) . iBagnulo et"al] d 19951) 
found the magnetic field in HR 1217 to be approximately dipolar, 
with a polar strength of Bp — 3.9 kG, an angle i = 137° between 
the line-of-sight and rotation axis, and angle /3 = 150° between 
the rotation and magnetic axes. More recent studies have confirmed 
that the field is primarily dipolar, with some variation in the value 
of the polar field strength, with [R yabchikova etal. (20^ finding 
Bp = 4.4 kG and lLuftinger et all ( l20IO() finding Bp = 4.2 kG. 

In this paper we seek to understand better how the amplitude 
of pulsation modes varies with time in HR 1217 through a wavelet 
analysis of the 1986 and 2000 data. Fu rthermore, since the analysis 
of the 1986 and 2000 data bv lKurtze t al. ( 1989) and Kurtz et al] 
d2005h did not use explicit statistical weighting, we also calculate 
weights for these data. Considering the wide range in the apertures, 
instrumentation, atmospheric and weather conditions, there is much 
to be gained by weighting the data appropriately. 



2 WEIGHTING SCHEME 

In weighting data, the estimated uncertainties, ai, are used. We can 
assign a weight, Wi, to the i"^ data point, as Wi = l/crf- The 
value of the exponent, x, depends on the nature of the noi se in the 
data: for Gaussian noise x — 2, which is commonly used. [Handle j 
([2003) found that a smaller exponent was more effective in reduc- 
ing the noise level for some data sets, so we do not assume x = 2, 
but rather find the value that minimizes the noise level. It must 
be noted that low-frequency sky-transparency noise was removed 
from the data by successive removal of low-frequency peaks (gen- 
erally below 0.6 mHz, and always well separated from the oscilla- 
tion signal at 2.7 mHz) to the point that low-frequency noise was at 
the same amplitude as noise at higher frequencies. This was done 
to obtain approximately white noise across the frequency spectrum 
so that least-square s errors could be ap prop riately estimated fo r the 
original analysis of lKurtz et all dl98Sd) and lKurtz et alj d2005l) . For 
our weighting scheme, this means that we do expect x ~ 2. 

For neither the 1986 nor the 2000 data do we have uncertain- 
ties available, and so we needed to determine them. To do this, the 
power in the range of frequencies in which oscillations are present 
was temporarily removed from the time series using a notch filter. 
The standard deviation, ai, in the residual time series was then used 
to calculate the weights. 

The timespan around each data point over which ai was cal- 
culated was treated as a free parameter. To find the best weighting 
scheme, we varied the timescale to find the one that minimized the 
average noise level at high frequencies with the exponent x = 2. 
Upon settling on this timescale, the exponent was varied to mini- 
mize the noise level further. 

The optimal timescale for both 1986 and 2000 was found to 
be about 20 minutes, a value presumably related to the typical 
timescale of variations in the atmosphere. The optimal exponent 
for the 2000 data was found to be 2.0, indicating that the noise dis- 



Variability in mode amplitudes in HR 1217 3 



15 

to 

E 5 
■S 

t-5 
E 

< 10 



^ 4 
to 

^ 2 



-I 



(a) 



; 1 



if! 



(b) 



; 1 



! 

i I 
I 



1 'I 



6720 6730 6740 6750 6760 

JD 2440000+ (days) 



6770 



6780 



Figure 1. (a) Time series of the photometric measurements of HR 1217 and (b) corresponding weights from 1986 October-December, colour coded as follows: 
South African Astronomical Observatory (SAAO, 1.0 m): red, (0.75 m): pink; European Southern Observatory (0.9 m): dark blue, (0.6 m): orange; Mount 
Stromlo and Siding Spring Observatory (0.6m): yellow; Lowell Observatory (1.1 m): black; Cerro Tololo InterAmerican Observatory (CTIO, 1.0m): dark 
green, (0.6m): light green; McDonald Observatory (McD, 2.1 m): maroon; Mount John Observatory (0.6m): cyan; Mauna Kea Observatory (MKO, 2.2m): 
purple, (0.6 m): magenta. The amplitude modulation with rotation is clearly visible in the top panel. 




Figure 2. Same as Fig.[T]but for observations from November-December 2000, colour coded as follows: CTIO (1.5 m): dark blue; SAAO (1.9 m): red, 
(0.75 m): pink; Observatorio del Teide de Instituto de Astrofisica de Canaria (0.8 m): black; MKO (0.6m): magenta; McD (2.1m): cyan; Perth Observatory 
(0.6 m): green; Beijing Astronomical Observatory (0.85 m): yellow; Aryabhatta Research Institute of Observational Sciences ( 1 .0 m): orange. 
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Figure 3. Power spectrum without weights {top) and with noise-optimized 
weights {bottom) between 1 and 9 mHz for the 2000 data. Clearly visible 
are the principal frequencies around 2.7 mHz. The use of weights has re- 
sulted in a reduction of noise of approximately 31% in amplitude allow- 
ing the second harmonic of the principal frequencies around 5.4 mHz to be 
distinguished above the noise level in the bottom panel. The peak around 
8.4 mHz, which at first may appear to be the third harmonic of some of the 
principal frequencies, is most likely instrumental in origin. 



9 mHz for the 2000 data, both with and without weights. It is clear 
that using weights has significantly improved the signal-to-noise, 
with the second harmonic frequencies around 5.4 mHz becoming 
distinguishable above the noise peaks. Also visible in this figure 
is a clear peak at 8.36 mHz. Interestingly, it falls at the third har- 
monic of the oscillations. (To be specific, it coincides with the third 
harmonic of the frequencies which had been missing in the 1986 
data.) However, the most likely explanation for this peak is a peri- 
odic drive eiTor as the telescopes track the star, given the period of 
2 sidereal min (corresponding to a frequency of 8.356 mHz). This 
8.4 mHz signal was more prominent at some telescopes, and not 
present at all at others, showing that it is instrumental in origin. 

A drawback of using weights is a worsening of the spectral 
window. This arises from an effective widening of the gaps between 
observations, which results from giving higher weight to superior 
data while down-weighting the beginnings and ends of nights, and 
from weighting small telescopes less compared to larger telescopes. 
In Fig.|4]the spectral window is shown for the 2000 data, with and 
without weights. The daily sidelobes have increased in amplitude 
by about 25 per cent. However, since these are multisite observa- 
tions, the daily sidelobes were already relatively small and the in- 
crease is not seriously detrimental: the reduction in the noise level 
more than compensates for the increased sidelobes when it comes 
to detecting weak peaks, as evinced by the second-harmonic fre- 
quencies in Fig.|3] 



3 RESULTS 
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Figure 4. Spectral window without weights {left) and with noise-optimized 
weights (right) for the 2000 campaign. 



tribution is close to Gaussian. For 1986 the exponent was slightly 
lower, at 1.8. The two time series and their optimal weights are 
shown in Figs[T]and[2] It is clear that the weights vary considerably 
over a night and between telescopes. 

With these weighting schemes, the noise at high frequencies 
in the amplitude spectrum was reduced by 22 per cent for the 1986 
data and 3 1 per cent for 2000. As an example of the improvement, 
Fig.[3] shows the noise level in the power spectrum between 1 and 



The amplitude spectra of the 1986 and 2000 data sets are compared 
in Fig.|5] The smoothed amplitude spectra (red lines) indicate the 
background noise level and accentuate the presence of significant 
peaks. Note the change in t he amplitudes of th e principal frequen- 
cies between 1986 and 2000. lKurtz etai] j2005l) cited this change as 
the reason that the frequencies vq and were not detected in 1986. 
In spite of our improvement in the noise level, these modes are still 
not clearly distinguishable above the noise level in the 1986 data. 
Mode frequencies were determined by iterative sine- wave fitting, 
yielding results in agreement with those reported by iKurtz et aU 
( 120051) . No additional modes were found. 



3.1 Wavelet Analysis 

To examine the variation in amplitude of the pulsation modes we 
performed a wavelet analysis. This involved calculating the Fourier 
transform of the time series inside a moving Gaussian window of 
fixed width (e.g. 20 h). We multiplied the weights of the data by 
this Gaussian window and then took the Fourier transform to reveal 
the strength of the peaks during this period. The Gaussian was then 
centred on a later time (e.g. 2 h later) and the process repeated. This 
resulted in a series of Fourier transforms that track the changes in 
amplitude of the peaks. 

By applying the Gaussian window to the weights, and not 
the observations, we maintained the calibration of the amplitudes. 
This is verified by the wavelet analysis of the spectral window 
(Fig.m bottom for the 2000 data). Although the amplitude of the 
central frequency remains constant, it is clear that there is a large 
amount of variation in the sidelobes due to the varying data cov- 
erage. This could potentially result in temporal variations in the 
amplitude of adjacent modes, especially since frequencies are ap- 
proximately separated by 3d~^ in HR1217. However, while the 
first and second sidelobes are heavily affected, the third sidelobe is 



Variability in mode amplitudes in HR 1217 5 




Frequency (mHz) 



Figure 5. Amplitude spectram of HR 1217 in 19 86 October-Decem ber (top) and 2000 November-December (bottom), unsmoothed (thin black line) and 
smoothed (thick red line). Frequencies as given bv lKurtz et all )2005h are marked by blue dashed lines (u\ to vg„ from left to right), with their first sidelobes 
marked by red dotted lines. 



barely affected, except towards the extremities of the time series. 
Keeping this in mind, we can reasonably expect that any variation 
shown at the frequencies of HR 1217 is due to real variation in the 
surface amplitude of the mode, resulting from temporal variations 
in either the intrinsic mode energy or in the form of the horizontal 
variations in the outer layers of the stellar envelope that cause the 
eigenfunctions to deviate from spherical-harmonic forms. 

Although windowing with a Gaussian reduces the effective 
length of the time series, and thereby permits more rapid changes 
to be observed, it causes a degradation of the frequency resolu- 
tion. The width of the Gaussian must be chosen appropriately to 
balance the time resolution against frequency resolution. Shorten- 
ing the time series has the added side-effect of increasing the noise 
level. 

Fig.[6](top) shows the results of the wavelet analysis on the 
data from 2000, using a Gaussian width of 20 h. Despite the fre- 
quent gaps in observations, which are indicated in this figure as 
white regions, the amplitude modulation of the principal frequen- 
cies is clearly apparent. This amplitude modulation arises because 
different aspects of the modes are visible as t he star rotates. The 
ephemeris prediction of amplitude maxi ma by Kurtz et al. 
based on the rotation period of 12.4572 d ( lKurtz&Marandll987h . 
is marked by the arrows. This ephemeris matches the observed 
maxima remarkably well, fourteen years after the observations 
on which i t was based. Principal frequencies, as determined by 
iKurtz et al. are marked by black horizontal lines. Other fre- 

quencies with significant power are aliases of the principal frequen- 
cies, which shift slightly in frequency as the sampling of the data 
within the Gaussian window changes. 



In the 1986 data, the gaps in the observations are longer and 
more frequent than in the 2000 data, with a comparable observation 
time, but spread over twice as long. The wavelet analysis for the 
1986 data was therefore done with a wider Gaussian (48 h). The 
false-colour plots are not shown here since it is difficult for the eye 
to follow changes across the numerous wide gaps. To look at the 
amplitude variation then in more detail we use an alternative way 
of displaying the results with a cross-sectional plot of the power 
contours, which shows the change in p eak power as a fun ction of 
time for each frequency determined bv lKurtz et al.l ( l2005h . This is 
shown in Fig.|7]for the principal frequencies in the 1986 data. The 
errors bars for each point were determined from the average noise 
level between 6 and 7 mHz. The power has been normalized in this 
figure, with the maximum power for each frequency over this time 
period being equal to 1. We do this to emphasize the change in 
strength of each individual mode over a rotation period. For clarity 
the plots have been offset vertically. 

The cross-sectional plot clearly shows the change in ampli- 
tudes over a rotation period. We also see that the rotational modula- 
tion of pulsation amplitudes is not completely regular. For the first 
pulsation maximum, all frequencies peak together. At the second 
maximum, v\ peaks before and after the other frequencies and is at 
a minimum when the other frequencies reach their maximum. This 
is similarly seen at the third maximum for v\ and vs- Note that these 
are the weakest frequencies. The epochs of maximum amplitude 
show evidence of a phase variation in the amplitude modulation, 
with one of the strongest frequencies, Vi, reaching maximum am- 
plitude later than the other frequencies on the third maximum. In- 
terestingly, these phase shifts are not the same at each maximum. A 
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Figure 6. Top Wavelet analysis showing power in the principal frequencies as a f unction of time for the 2000 data. White regions indicate an absence of 
observations. Black horizontal lines i ndicate the p rincipal frequencies as given bv lKurtzetai] )2005h and aiTows indicate the ephemeris prediction of the 
maximum of pulsation amplitude from lKurtz et al.1 fl989i) . The power scale has been truncated to show lower amplitude frequencies. Bottom Contour plot of 
the spectral window of the 2000 data. Black horizontal lines indicate the approximate mode spacing ( Ai/ = 34 /^Hz). 



third observation is the considerable variation in the power of each 
frequency at their maxima, with 1^2 showing the greatest change be- 
tween the power observed at the first maximum, and at the second 
and third. 



The cross-section of the 2000 data shows further evidence 
of these behaviours. The normalized cross-section of the princi- 
pal frequencies is shown in Fig.[8] The ephemeris prediction of 
iKurtz etal.ll 19891) is shown by vertical lines. Once again, the am- 
plitude modulation with rotation is iiTegular. Many of the frequen- 
cies show multiple peaks at one or both maxima (for example 1^5 
and vn), or plateau instead of having a distinct peak (in particular 
U2). No frequency shows the same structure at both maxima. The 
phase difference of the maxima between different modes is most 
clearly apparent at the second maximum between uz, U4 and v^,. 
The variation in the power of each mode from cycle to cycle that 
is also visi ble is most likely the source of additional frequencies 
reported bv lKurtzetai] ( l2005h that were close to the principal fre- 
quencies, and not separated by the rotation frequency. 



4 DISCUSSION 

The wavelet analysis of HR 1217 has revealed interesting behaviour 
in the form of structure in the rotational modulation, phase differ- 
ences in the rotational modulation between modes and variation of 
the pulsational energy in each mode. We consider here possible ex- 
planations for these observations. 

Several theoretical studies investigating the interaction 
between the magnetic field and pulsations in roAp stars have 
found that the pulsation modes are severely distorted from a 
single spherically harmonic structure jpz iembowski & Goo(3 
Il99d ICunha & Goughl I2OO0I : ISaio & Gaiitschvll2004b . Distortion 
can also result from lateral variation of the vertical stratifica- 
tion of the star that o n e expe cts to be associated with spots 
I Montgomery & Gougri l2003h . i aio. Rvabchikova & Sachkovl 
1 20100 modelled the pulsations of HR1217, with the oscillation 
frequencies identified as alternating between modified 1=1 and 
^" = 2 modes, with the irregularly spaced vg, identified as a modified 
£ = 3 mode. They have labelled these modified modes as Im = 0, 
1, 2 and 3, corresponding to the I value of the component that has 
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Figure 7. Power of principal frequencies as a function of time for the 1986 data, showing the amplitude modulation with the rotation period, normalized such 
that the peak power duiing the two central amplitude maxima is equal to one. Error bars indicate the average noise level between 6 and 7 mHz. The dashed 
lines c onnect the data points to 'guide the eye'. The black vertical lines are the ephemeris prediction of the maximum of pulsation amplitude from lKurtz et alj 
<1989l) . 
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Figure 8. Same as Fig.|2] but for the 2000 data, the peak power during the two central amplitude maxima is normalized to one. 
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the greatest kinetic energy in the expansion of the eigenfunctions 
using the axisymmetric spherical harmonics Yf™^". 

The amplitude modulation with rotation phase for modes 
of Im = 0, 1, 2 and 3 has be en predicted (see Fig. 8 of 
ISaio. Rvabchikova & SachkovlbOlOh . They assumed that the pul- 
sation and magnetic axes in HR 1217 are aligned, with the mag- 
netic geometr y given by the param eters i = 137°, /? = 150°, as 
determined bv lSagnulo et al.l ( ll995[) . A single peak at the pulsation 
maximum, in phase with the magnetic maximum, is expected for 
the modified dipole, quadmpole and octopole modes. For modified 
radial modes the amplitude modulation is smaller, with the pulsa- 
tion maximum expected to occur at magnetic minimum. Does this 
theoretical model well explain the structure in the observed ampli- 
tude modulation? 

All observed frequencies exhibit an amplitude modulation 
which, in general, peaks around the magnetic maximum, consis- 
tent with their being modes of /,„ = 1, 2 or 3. However, more 
detailed structure is also observed. The modulation of v\ and 
at some rotational cycles in 1986 is particularly hard to explain: a 
pulsation minimum is observed at magnetic maximum. While this 
might seem to imply a radial mode geometry, pulsation maximum 
is not observed at magnetic minimum, but rather at the first and 
third quarter of the magnetic cycle, which is inconsistent with a 
lm = Q mode. 

The best-fit mo del of HR 1217 by 

ISaio, Rvabchikova & Sachkovl ( 1201 Ol) does have = and 
Im = ^ modes closely spaced in frequency, and v\ is identified 
as an = 1 mode in this model. It is conceivable that v\ might 
behave as a combination of lm=Q and /,„ = 1 modes. However, 
a linear combination of the Im = Q and Im = 1 modulation does 
not suffice to explain the observed modulation for u\ and z/g. 
Furthermore, vs, is identified as an = 3 mode in this model, 
and is well separated in frequency from the nearest radial mode. 
Observed structure in the modulation that is more complex is 
even more difficult to explain by mode geometry distortion. If this 
is indeed the correct explanation, then it would appear that the 
distortion to the mode geometry can change over a remarkably 
short period of time, and a new mechanism would be required to 
explain this. 

It is already well known that there is a phase difference be- 
tween radial velocity variations and luminosity variations, based 
on simult aneous spectroscopic and photometric obser vations of 
HR 1217 dSachkov et alj|2004 [Rvabchikova et alj|2007l) an d other 
roAp stars (e.g.l Sachkov et al.l2 008VMkrt ichian et al.l2008h . There 
are also phase differences in the radial velocity variations of absorp- 
tion lines of different elements which, due to chemical stratification 
in the atmosphere of roAp stars, may be attributed to depth effects 
jRvabchikova et aLll2007h . However, the rotational phase lag ob- 
served between different modes in the wavelet analysis is different: 
a lag in the rotational modulation as opposed to a lag in the pulsa- 
tions themselves. 

Magnetic Doppler imaging of HR 1217 by iLilftinger et al.l 

( I2OIOI) found quite interesting elemental abundance patterns. Re- 
gions of abundance enhancement or depletion were seen around 
either the rotation phase when the positive magnetic pole is visible, 
or where the magnetic equatorial region dominates the visible sur- 
face, depending on the chemical element. Most interesting is that 
the enhancements for different elements were shifted in longitude 
relative to each other. 

Since the observations in 1986 and 2000 were in the photo- 
metric B band, a range of atmospheric depths was being observed. 
The different modes have slightly different pulsation cavities, and 
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Figure 9. The total power in all modes observed in 1986 {vi — , ug) (top) 
and 2000 (ui - us) (bottom). 



amplitude and phase variations are strongly depth dependent in 
roAp stars in general, so it is plausible that this could explain the 
phase behaviour. Once again, it seems that the phase difference can 
change over a short period of time. 

It is apparent from a comparison of the amplitudes at maxi- 
mum in both Figs|7]and[8]that, in addition to the variation in power 
due to rotational modulation, the intrinsic amplitude of each mode 
varies over as short a timescale as one rotational period (12.4572 d). 
Short-timescale va riations have been observed in oth er roAp stars, 
such as HD 60435 l lMatthews. Wehlau & KurtJ 19871) . It is difficult 
to resolve power variations on even shorter timescales for HR 1217 
owing to the modulation caused by rotation. Could these power 
variations on timescales shorter than the rotation period be respon- 
sible for the behaviours that have been observed, namely structure 
in the rotation modulation and phase difference between modes? 
If the pulsation energy in a mode varies considerably over a few 
days then it is conceivable that this could result in the behaviour 
seen. We are then left, however, with several coincidences, such as 
and vs being at a minimum at both full and half rotation phase 
occasionally in 1986, and a feature in 2000 in which, at the first 
pulsation maxima, several frequencies appear to show two distinct 
maxima. While a change of power in the modes is certainly re- 
sponsible for some of the observed features, it seems likely there 

are other processes involved as well. 

As previously mentioned, iKurtz et al] ( l2005h speculated that 
the total pulsational energy could be conserved, with nonlinear in- 
teractions transferring energy between modes. Fig.|9]shows the to- 
tal power from all modes as a function of time. While it can be 
seen that the power at maximum reaches approximately the same 
value at several maxima, both in 1986 and 2000, there are two no- 
table exceptions, specifically the third maximum in 1986 which 
only reaches just over half the height of the other maxima, and the 
first maximum in 2000, which exhibits the previously mentioned 
double peak. Caution is required in relating pulsation power to the 
pulsation energy, the determination of which requires a complete 
characterisation of the modes. 
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5 CONCLUSIONS 

We have reanalysed the 1986 and 2000 multisite photometric obser- 
vations of theroAp starHR 1217, using a weighting scheme to min- 
imize the noise level. Despite the improvement in the noise level, 
the 'missing' frequencies i^e and U7 are still not detectable in the 
1986 data, and no other unreported frequencies were discovered in 
either data set. 

A wavelet analysis of the 1986 and 2000 data has revealed the 
rotational modulation and temporal variation of individual modes. 
The rotational modulation of all modes is, in general, consistent 
with the exp ected modulation for modified 1, 2 and 3 modes, as 
predicted bv lSaio. Rvabchikova & Sachkovl ( I2OI0I) . However, this 
leaves the problem of the small spacing, Si>i2 = v{n, Zm = 1) — 
0.5[i'{n, Im = 2) + fin — l,lm = 2)], which, if the Im = 1 and 
= 2 identifications are correct, is as small as 0.5 /iHz. For a ll 
of the models calculated bv lSaio, Rvabchikova & Sachkovl l l2oTol) . 
this spacing is always larger than ~ 3 jjHz. A solution to this, in 
which the = 1 modes are, in fact, Im ~ modes, is ruled out 
by the rotational modulation observed in this wavelet analysis. 

Temporal variations of modes are also seen in the wavelet 
analysis. These fall into three broad categories: variation in the pro- 
file of the rotational modulation, variation in the phase of maximal 
power and variation of the power itself. The explanation of these 
observations remains an open question. We have considered sev- 
eral possibilities. A variation of the geometry of some modes may 
explain variations in the profile, although it does seem unlikely that 
this could fully explain the variations seen. Depth effects could re- 
sult in the phase variations observed between different modes. All 
of the behaviours could be explained, at least in part, by a variation 
in the pulsational energy of each mode, although a cogent theory 
to explain the deta ils of such variation is lacking. As suggested by 
iKurtz et al. I j2005l) . there may be some mechanism whereby some 
of the energy of one mode may be transferred to another. Of par- 
ticular note for all of these behaviours is the short timescale over 
which these variations are seen to occur. 

Despite the high quality of the 1986 and 2000 data, improved 
further by our use of weights, multisite ground-based campaigns 
are still limited by gaps in the data. Observations of HR 1217 from 
a dedicated asteroseismic satellite would be useful. Such obser- 
vations do exist: the Mircmvariability and Oscillations of Stars 
(MOST) satellite observed the optical intensity of HR 1217 for 
thirty days in November-December 2004. The MOST observations 
are taken in white light and are thus not directly comparable to 
these ground based observations taken in the B band since they 
sample a wider range of atmospheric depth. Nevertheless, a simi- 
lar analysis would prove enlightening. Much more may be gained 
from spectroscopic observations. The chemical stratification of the 
atmosphere of roAp stars reveals the radial structure and level of 
nodes of the modes and their harmonics, revealing the effect of at- 
mospheric depth in this analysis. Our results show that much of the 
physics of roAp stars is still not well understood. More advanced 
models are required that can reproduce variations in the pulsational 
energy of different modes. 
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